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ABSTRACT 

The conservation and convergence properties of spectral Fourier methods for the numerical 
approximation of the Korteweg-de Vries equation are analyzed. It is proven that the (aliased) 
collocation-pseudospectral method enjoys the same convergence properties as the spectral Galerkin 
method, which is less effective from the computational point of view. This result provides a precise 
mathematical answer to a question raised by several authors in the latest years. 
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In this paper, we analyze the numerical approximation by Fourier spectral methods to the 
Korteweg-de Vries ( briefly K.d.V. ) equation with periodic solutions: 
au/at + u du/dx + cx d 3 u/dx 3 = 0 , x € IR , t>0 , 

(1.1) u(x+2n,t) = u(x,t) , x € [R , t > 0 , 

u(x,0) = u°(x) , x € IR , 

where the initial condition u° is a real valued 2n-periodic function and a is a real non zero 
parameter. 

The Korteweg-de Vries equation was formerly introduced in 1 895 by Korteweg and de Vries 
[10] to model long, unidirectional, dispersive waves and, nowadays, it is used to describe 
phenomena arising from several different fields of applied science. A survey of these numerous 
applications is given, e.g. , in Bardos [ 1 ]. 

From a more theoretical point of view, after the early sixties, the study of the K.d.V. 
equation has benefitted from the application of the scattering theory and from the discovery of 
many energy integrals. We refer to the pioneering work by Miura [14], Temam [22], and the 
more recent papers of Lax [9] , Miura [ 1 5] and Bardos [ 1 ,2], Moreover, the analogies with the 
studies of Hamiltonian systems that have been focussed lately, have permitted extension of the 
applicability of the K.d.V. equation to new theoretical and physical models. These aspects, and their 
relevance for the interaction between pure and applied mathematics, are discussed in a very 
instructive review article of Bardos [ 1 ]. 

Numerical approximations of the K.d.V. equation based either on finite differences or finite 
elements methods are abundant in the literature. We refer, interalia, to the papers by Bona, 
Pritchard and Scott[4] , Bona, Dougalis and Karakashian [3] and to the references quoted therein. 
Fourier spectral methods have been used also in many applications of the K.d.V. equations in the 
last decade. We refer, e.g., to the works by Tappert [21], Fornberg [5], Schamel and Elsasser 
[20] , Fornberg and Whitham [6] , He Ping and Ben Yu [ 1 3] and Pasciak [18]. 

The classical Fourier-Galerkin method has been used, as well as the more flexible 
collocation pseudospectral method in which the discrete Fourier transform is applied to deal with 
nonlinear terms ( see e.g. Canuto, Hussaini, Quarteroni and Zang [5], chapters 2 and 4 ). A crucial 
question has been raised from several authors, and this is precisely whether the (aliased) 
collocation-pseudospectral method retains the same asymptotic accuracy as the pure Galerkin 
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method. This is a sort of master question which is recurrent in the context of numerical 
approximations by the spectral method. In the case of the K.d.V. equation, particularly, this doubt 
has induced several authors to introduce new pseudospectral methods with the aim of regaining the 
(presumably lost) exponential accuracy ( see again Canuto, Hussaini, Quarteroni and Zang [5], 
sections 4.5 and 4.6 ).ln this paper, we provide a precise mathematical answer to the above 
question. We prove that the genuine (non dealiased, non skew-symmetric) collocation 
-pseudospectral method enjoys the same convergence properties as the Galerkin method. 

In section II we start by proving that the Galerkin approximation conserves the three first 
energy integrals of the K.d.V. equation. Then, classical energy methods allow us to prove that the 
Galerkin solution converges with spectral accuracy to the mathematical solution. 

In section III, the collocation-pseudospectral approximation is considered. The related 
solution fails to conserve the second energy integral. However , convergence with spectral accuracy 
in any finite time interval [0,T] can still be shown by using a much more involved proof. In its 
essence, the new proof exhibits first that the pseudospectral solution cannot blow-up in a small 
subinterval [O.t,] of [O.T], This property, joined to the property of consistence of the numerical 
method, allows us to initialize an inductive process which yields the desired result on the large 
time interval [O.T], 

In this paper, we will not be concerned with any time discretization of the K.d.V. equation. 
However, we recall that the semi-implicit time advancing schemes are customarily used for such 
a kind of equation. These schemes are computationally convenient since, at each time interval , they 
yield a diagonal system in terms of the unknown Fourier coefficients of the spectral solution. 
Moreover, for finite time intervals, they are stable without any restriction on the time and space 
discretization parameters. We refer the interested reader to Chan and Kerkhoven [6] where a 
linear stability analysis is presented for the K.d.V. equation, to Quarteroni [ 1 9] where a nonlinear 
stability analysis is carried out for a family of equations of the same kind and to Bona, Dougalis and 
Karakashian [3] and the references quoted therein where extended equations are also considered. 

Working along with 2n-periodic functions, we introduce the periodic Sobolev spaces 
defined over ]0,2n[ . We first recall the definition of classical Sobolev spaces. We set 
L 2 (0,2tt) = { f: ]0,2n[-*-C , II f II = [\l* |f(x)| 2 dx ] 1/2 < oo } , 
and we denote its scalar product by (.,.). Now , for any integer r > 0 , we set 
H r (0,2n) = {f€ L 2 (0,2n) , II f || r = [Lj \ 0 l|a j f/dx j || 2 ] 1/2 < oo } , 
and for any real r > 0, not in DM, the space H r (0,2n) is defined by interpolation between 
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H E<r) (0,2n) and H E(r,+ 1 (0,2n) ( we denote here by E(r) the integral part of r ). Next we consider 
the subspace C"(0,2n) of C°°(0,2n) of all functions that are 2n-periodic as well as all their 
derivatives. Moreover, for any real r > 0, H^(0,2n) stands then for the closure of C“(0,2n) in 
H r (0,2n). As pointed out in Lions and Magenes [12], if r— 1/2 is not an integer, then H^(0,2n) 
consists of all functions of H r (0,2n) that are 2n-periodic as well as their derivatives of order 
^ r- 1 /2 . For any real number r ^ 0, we define H r (Q) as the dual space of H" r (Q). Its norm is 
again denoted by || . || r . 

Finally, if A is an interval of (R and X is a Banach space, for any function f from 1R into X, we will 
define 

II f Hl°° (A X) = SU P II «0 lx 
t € A 

With these notations it is known (see Temam [22]) that as soon as u° belongs to H™(0,2 tt), with 
m in DM, then the solution to (1.1) satisfies 

(1.2) II u IIl°° (o,t ;H™(0,2tt)) ^ ^ U ° ’ 

for any T > 0, where the constants are independent of u°. 

It is well known, that the family 

(1.3) <p k (x) = (2n) -1/2 exp (ikx) , k € Z , 

is orthonormal and complete in L 2 (0,2n). Thus a natural approximation of L 2 (0,2n) by periodic 
functions will consist of the spaces defined by 

(1.4) V N € 2(N , S N = span { <p k , -N/2 ^ k < N/2 } . 

Let us denote by P N the operator 

o _N/2 

V g € L (0 ,2n) , P N g = 2Ik=- n /2 9k *Pk ■ 

with 

(1.5) g k = Jo n g(x)^p k (x) dx , k€Z . 

Since P N is in fact the orthogonal projection operator over S N we have equivalently 

(1.6) VimS h , fo n (P N g-g) «l> dx = 0 

For all g in L 2 (0,2n), (P N g) NeW converges to g. Moreover, it can be proved that for any r ^ s , 
r ^ 0 one has ( see, e g. Jackson [9], Pasciak [17]) 

(1.7) V g € H^(0,2n) , || g - P N g || s ^ cN s ~ r || g || r . 

(Throughout this paper, c will denote a positive constant, independent of N, not necessarily the 
same in different contexts). 
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Besides, all the norms defined bg the imbedding of S N into hf (0,2n) are equivalent since $ N is a 
finite dimensional subspace of H^(0,2n); more precisely we can readily see that 

(1.8) V (r,s) € (IR + ) 2 with r ^ s , V <p € S N , || <p || r < || <p || 8 < *(s) N s ' r || <p || r . 

The second inequality is known as the "inverse inequality". 

We finally notice that 

(1.9) V g € S N , || g || 2 = L^n /2 ( 1 + k 2 ) 9^ < 9o + 2 II dg/dx || 2 
This property will be frequently used in the sequel. 


A spatial approximation (continuous in time) of problem (1.1) based on the 
Fourier-Galerkin method reads as follows : 

Find a mapping u N : [O.T] — ► S N such that 

(11.1) V € S N , Vt, 0 < t < T , (9u N /9t + u N 0u N /0x + a 9 3 u N /9x 3 , hi) = 0 , 
u N (0) = P N u° . 

This entails a nonlinear system of O.D.E.'s for the Fourier coefficients (u N ) k (t) of the solution u N . 
We present now the main properties satisfied by the above Fourier-Galerkin approximation. They 
are concerned with the concepts of conservation, stability, uniqueness and convergence. 

Lemma 11.1 : There exists a unique solution u N to problem (11.1). Moreover this solution 
conserves the three first energy integrals of the K.d.V. equation, namely 

(11.2) (9/9t) [Jo" u N (x,t) dx] = 0 , 

(11.3) (9/9t) [Jo" lu N (x,t)| 2 dx] = 0 , 

(11.4) (9/90 [Jo" (a(0u N /9x) 2 - u 3 /3)(x,t) dx] = 0 . 

Proof : The existence of a maximal time t 0 , 0 < t 0 ^ T such that, for all t < t 0 , there exists a 
unique solution u N (t) to problem (11.1) is a classical result of the theory of differential systems. 
The problem is to get the existence for a "long" time T, or equivalently to prove that one can take 
t 0 = T. This result will be achieved with the use of (11.3) that the solution cannot blow-up. 

Since the initial condition u° is real, then u N (0) is real too. We deduce that u N (t) is real for 
any t < t 0 from the uniqueness of the solution to problem (11.1). 

In order to show (11.2), let us first choose qi = 1 as a test function in (11.1). We get 

(9/9t)Jo* u N (x,t) dx + ( 1 / 2 ) 10 * ( au N /ax ) ( x .O dx + <* Jo" (d 3 u N /9x 3 ) (x.t) dx = 0 , 
using the periodicity of u N , we deduce then (11.2). Choosing now hi = u N in (II. 1 ) , we obtain 

(11.5) (d/dOjo* (u N (x.t)) 2 dx +( 1/3)Jo* (du^/dx) (x,t) dx + ocJq* (u N 9 3 u N /9x 3 ) (x,t) dx = 0 
Integrating by parts and using the periodicity of u N yield 

Jo* (u N 9 3 u n /9x 3 ) (x ,t) dx = - ( 1 /2) Jo* 0/9x (9u N /9x) 2 (x ,t) dx = 0 . 

Similarly we have 

Jo"(9u 3 /9x)(x,t)dx = 0 . 

We derive now (11.3) from (11.5). Integrating (11.3) between 0 and t 0 proves that no blow-up 
occurs at time t 0 . More precisely for any t, 0 < t < t 0 , we derive 

(11.6) II U N ( . ,t) II D II U„( . ,0) II < II U° II 
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whence t 0 is equal to T and we can state that the existence and uniqueness of the solution u N holds 
for any time t, 0 < t < T . 

In order to prove now (11.4), let us take 41 = P N [u^ + 2a 9 2 u N /3x 2 ](. ,t) in (11.1) and, for 
convenience of notation, let us drop for a while the explicit dependence on x and t. Then 

(11.7) Jo” (Ou N /Ot) P N [u 2 + 2a 9 2 u N /9x 2 ] dx 

+ ( 1 /2) Jq 11 (9/9x)[u^ + 2a d 2 u N /dx 2 ] P N [u^ + 2a d 2 u N /dx 2 ] dx = 0 . 

Since du N /dt is in S N we have from (1.6) 

Jo” (au N /at) P N [u^ + 2a 9 2 u N /9x 2 ] dx = Jo" 9u N /at [u 2 + 2a a 2 u N /3x 2 ] dx 

=(a/at) [Jo" (u^/3 - a(au N /ax) 2 ) dx] . 

On the other hand using again (1.6) we derive 

Jo”(9/9x)[u^ + 2aa 2 u N /ax 2 ]P N [u^+2aa 2 u N /ax 2 ]dx 

= ( 1 /2)Jo"(9/9x) {P N [u 2 +2a9 2 u N /9x 2 ]} 2 dx □ 0 . 

From (11.7) we get now (11.4 ). 

Remark 11.1 : The estimates (11.2) to (11.4) are the discrete analogues of the conservation laws 
for the K.d.V. equation defined in, e.g. , Miura, Gardner and Kruskal [16] which are at the basis of 
the scattering theory. 

In the next two lemmas we state some a priori estimates for the Fourier-Galerkin solution 
in higher order norms. 

Lemma 11.2 : Assume that u° belongs to (0,2n). Then there exists a constant c > 0 independent 
of N such that for any t , 0 < t < T: 

(11.8) || u N (.,t) |,< c. 

Proof : For any t < T, we derive from (114) 

(11.9) Jo” (a(9u N /9x) 2 - u^/3)(x,t) dx = Jq" (a(9u N /9x) 2 - u^/3)(x,0) dx . 

Using now the continuous imbedding of hHo^tt) into L°°(0 ,2tx) ( see, e.g. [ 1 2]) , we obtain first 
that 

Jo”(uj|/3)(x,0) dx <(1/3) || u N (.,0) || L oo || u N (. ,0) || 2 
<c| u N (.,0) ||, || u N (.,0)|| 2 
By the definition of u N (. ,0) and (1.7) it follows that 

(11.10) fj”(u 3 /3)(x,0)dx < c II u° II, ||u° II 2 < c ( 11 u° Ilf * ||u° ||“ ) 

In a similar way, using (11.6), we get 
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Jo" (ujj/3)(x,t) dx < o 1 u N (.,t) II, n u N (.,t) II 2 *S c II u„(.,t) II, II u° f , 

<<W/2)|u,(.,OI? * o » u° || 4 , 

The estimate (11.8) is then an easy consequence of (1.9), (11.9) and (II. 1 0) . 

With this stability in the H^(0,2n)-norm we can prove now, as in the continuous case, the 
boundedness of another energy integral. 

Lemma 11.3 : Assume that u° belongs to H 2 (0,2tt). Then there exists a constant c > 0 independent 
of N such that for ang t, 0 < t < T: 

( 11 . 11 ) ||u N (.,t)|| 2 <c. 

Proof : As in the proof of Lemma 11.1 we have to choose properly a test function in (11 1 ). This 
time we take vp = P N [u^ + 3oc(9u N /9x ) 2 + 6ocu N 9 2 u N /9x 2 + (18/5)oc 2 9 4 u N /9x 4 ](. ,t). This 
choice yields ( here again we drop the dependence on x and t ) 

(11.1 2) Jo" [9u N /9t 4 < + u N (9u n /9x) y + a (9 3 u N /9x 3 ) 5j»] dx = 0. 

Let us examine now the first term in (11.1 2). 

Since 9u N /9t belongs to S N , by (1.6) we have 

{o^dUfj/dt dx = }Q n 9u N /9t(u^ + 3a(9u N /9x) 2 + 6au N 9 2 u N /9x 2 + ( 1 8/5)oc 2 9 4 u N /9x 4 )dx; 
integrations by parts yield 

(11.13) (o"9u N /9t 4 . dx = (9/9t) {*" [uJ/4 - 3cxu N (9u N /9x ) 2 + (9/5) a 2 (9 2 u N /9x 2 ) 2 ] dx . 

Now we notice that 4 i = P n [u^ + 3oc(9u n /9x) 2 + 6ocu n 9 2 u n /9x 2 ]+( 1 8/5) oc 2 9 4 u N /9x 4 , hence we 
obtain 

(II. 1 4) Jq"[ u n (9u n /9x) + oc9 3 u N /9x 3 ]vp dx=A + B+ C + D + E 
with 

A := ( 1 8/5)a 3 Jo" 9 3 u n /9x 3 9 4 u n /9x 4 dx = 0 , 

B := J 0 u N (9u N /9x) P N (uJ) dx , 

C:=aJ 0 2 "[9 3 u N /9x 3 u 3 + 3 u N 9u N /9x P N (u N 9 2 u N /9x 2 )] dx , 

D := 3ocJo" u N 9u n /9x P n ((9u n /9x) 2 + u N 9 2 u N /9x 2 ) dx , 

E := oc 2 Jo" [( 1 8/5) u N 9u n /9x 9 4 u n /9x 4 + 3 9 3 u N /9x 3 ((9u N /9x) 2 +2 u N 9 2 u N /9x 2 )] dx 
In order to bound B we use the continuous imbedding of H^(0,2n) into L°°(0,2n) 

B c ||u N ||, ||9u n /9x ||o || P N (u 3 ) || c ||u N ||? || u 3 || < c ||u N ||? || u 3 ||, 

Since H ],((), 2n) is an algebra, || uj ||, < c || u N || 3 and therefore we deduce from (11.8) that 
(11.15) | B | ^ c 

Let us consider now the term C; integrating by parts the first addendum we obtain 
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C = - a}^ [0 2 u n /3x 2 (3 u 2 au N /3x) -3u N 3u N /3x P N (u N 3 2 u N /3x 2 )] dx 
= 3ocjo u N 3 u n /3x [P N (u N 3 2 u N /ax 2 ) - (u n 3 2 u n /3x 2 )] dx , 

so that 

C < 3a || u N | L oo || u N ||, || P N (u N 3 2 u N /3x 2 ) - (u N 3 2 u N /3x 2 ) || 

From (1.7) ( using again the imbedding of H^(0,2 tt) into L°°(0,2n)), we derive 

I C U c || u N || L co || u N ||, || u n 3 2 u n /8x 2 || < c || u N || L oo || u N ||, || u N || L oo || 3 2 u N /3x 2 || 

< c || u N || 3 || u N || 2 , 
and by (118) we conclude that 
(11.16) | C | < c || u N || 2 . 

We now have 

D = 3ajQ n u N 3 u n /3x P N ((a/3x)( u N au N /3x)) dx , 

= 3oc|q P n (u n 3u N /dx) P n ((3/8x)( u n 3u N /ax)) dx 
= 3a|Q n P N (u N 3 u n /3x) (3/3x) P n ( u n 3u N /ax) dx = 0 . 

Similarly we obtain . 

E = cx 2 \l " [(- 1 8/5) ((8u n /3x) 2 + u N 8 2 u N /3x 2 )(3 3 u N /8x 3 ) 

+ 3 (3 3 u N /3x 3 )((8u N /ax) 2 +2 u n 3 2 u n /8x 2 )] dx 
= -a 2 [(3/5) (8u N /ax) 2 a 3 u N /8x 3 - (6/5) u N 3/3x(a 2 u N /ax 2 ) 2 ] dx 

= a 2 j 0 [(3/5) (a/8x)((3u N /ax) 2 )3 2 u N /3x 2 + (6/5) u N (a/3x)(a 2 u N /3x 2 ) 2 ] dx . 

After integration by parts, we deduce that E = 0 . 

From (11.14), (11.15), (11.16) and the fact that A = D = E = 0we derive 
| }q"[( 1 /2) QuJ/ax + oc 3 3 u n /3x 3 ]hi dx | ^ c ( || u n || 2 + 1 ) . 

Due to (II . 1 2) and (11.13) we obtain 

(3/at) Jo" [ujj /4 - 3au N (3u N /3x) 2 + (9/5) a 2 (3 2 u N /ax 2 ) 2 ] dx ^ c ( || u N || 2 + 1 ) . 

After integrating between 0 and t, 0 < t < T and using (11.8) to bound the two first resulting terms 
under the integral, we have 

Jq" (9/5) a 2 (3 2 u N (. ,t)/8x 2 ) 2 dx < c ( 1 + J* || u N (,,s)|| 2 ds) 
where c depends on the H 2 (0,2n)-norm of the initial condition u° and on T. Using now the 
Gronwall's lemma yields (II. 1 1 ). 

We turn now to the convergence estimate for the Galerkin approximation of the K.d.V. 


equation. 
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Theorem 11.1 : Assume that u° belongs to H™(0,2n) , for some integer m $s 2 . Then there exists 
a constant c > 0 independent of N such that for any t, 0 ^ t < T: 

(11.17) || u(..t) -u M (.,t)D<cN , - m . 

Proof : For any time t, 0 < t < T, we set e(t) = P N u(t) - u N (t). From (1.1) and (II. 1 ), and setting 
E[f,g] = fdf/dx - gdg/9x, we derive for any y in S N 

(11.18) (9e/9t + oc9 3 e/9x 3 , »v) = (E[P N u,u] - E[P N u, u N ] , y) . 

Let us choose = e as test function in (11.18) and bound each term on the right-hand side. We 
obtain first, using (1.2) and (1.7) 

(11.19) | (E[P n u,u] , e) | < (|| u (I, + || P N u ||, ) || u - P N u||, ||e|| < c|| u - P N u ||, ||e|| . 
Moreover, a repeated use of integration by parts yields 

|(E[P n u,u n ] , e) | =|(1/2) Jo"(9/9x)[(P n u) 2 - u N 2 ]edx| 

= |( 1 /2) j 0 (P N u + u N ) e 9e/9x dx | 

= K 1 /4) \ 2 q (9/9x)[P n u + u N ] e 2 dx | . 

From (1.2) , (II. 1 1 ) and the imbedding of H^(0,2n) into L°°(0,2n) we deduce that 

(11.20) | (E[P n u,u n ] , e) | ^ c ||(9/9x)[P n u + u N ] || L «|| e || 2 ^ c || e || 2 . 

Let us note now that e(0) = 0; by'( 1 1 . 1 8) to (11.20) and the Gronwall's lemma we obtain 
|| e(t) || c exp(c t) [ {q || u(s) - P N u(s) || 2 ds] 1/2 . 

The estimate (11.17) is then an easy consequence of (1.2) , (1.7) , and the triangle inequality: 

II u - u N || < II u - P N u II + II e|| . 

The above result yields the following error estimate in the iV (0,2n)-norm. 

Corollaru 11.1 : Assume that u° belongs to H™(0,2ti) , for some m > 2 .Then there exists a 
constant c > 0 independent of N such that for any t, 0 < t < T: 

(11.21) || u(. ,t) - u N (. ,t) ||, < c N 2_m . 

Proof : This result is a consequence of the inverse inequality in (1.8). Indeed, we have 
II u N ( ,t) - P N u(.,t) I, <Nflu M (.,t)-P N u(..t)|| 

< N [ H u(.,t) - u N (.,t) || + || u(.,t) - P N u(. ,t) || ] , 
and from (1.7) and (11.17) we derive 
|| u N (,,t) - P N u(.,t) ||, ^ N 2_m . 

Now (11.21) follows by using again (1.7). 
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III. ANALYSIS OF THE COLLOCATION METHOD FOR THE APPROXIMATION OF THE K.d.V 

EQU AI lQkL 

Despite its mathematical interest, the Fourier-Galerkin method is generally abandoned in 
the applications in favor of the Fourier collocation method. The latter method allows a very 
efficient treatment of the nonlinear term u(du/dx) by transform techniques at the expense, 
however, of introducing an extra error due to the aliasing. This has induced many authors to 
dealiase the Fourier collocation solution by resorting to different kinds of techniques. For a more 
involved discussion about these arguments, the reader is referred to [5] ( see e.g. Chapter 3 and 
section 4.4.2 ). 

We will show in this chapter that the aliased Fourier collocation method is stable and 
convergent, and that its asymptotic rate of convergence is the same as that of the Fourier Galerkin 
method. 

111.1 Position of the problem. 

Let us introduce the collocation points Cj = 2nj/(N+1), for ] = 0.....N. Then we associate 
with this set the interpolation operator i N : C°(0,2n) -► S N defined by: 

(111.1) V f € C°(0,2n) , l N f€S N andI N f(Cj) =f(Cj) , j = 0,...,N. 

Now we define the pseudo-spectral derivative operator a N as 

(in. 2) v f € cVo^n) , a N f = a o i N f = a(i N f)/ax. 

Remark 111.1 : The calculation of the nodal values of a N f in terms of those of f is accomplished by 
two F.F.T.s plus N complex multiplications ( see, e.g., [5] Chapter 2 ). This requires 5N log 2 N 
operations if N is a power of 2. 

Let us introduce the following "discrete" scalar product over C (0 ,2n) 

(111. 3) V ip , ip € C°(0 ,2n) , («p,» t .) N = ((2n)/(N + 1))L J N =0 «p(t; j )^C j ) . 

It is well-known that it coincides with the L 2 -scalar product when the product <pip belongs to S 2N , 
hence in particular 

(111. 4) V <p , »p € S N , (ip ,»p) N = (<p ,ip) . 

Then the operator l N is precisely the orthogonal projection operator onto S N with respect to 
(.,.) N . Moreover.it has been proved in [17] that this operator satisfies the following inequality 
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for any real numbers r and s , 0 < r ^ s , 

(111. 5) V <p € H’(0,2n) , || <p - I N <p || r < c N r - S || «p || s . 

With these notations we can now introduce the formulation of the approximate problem 
obtained by a collocation pseudo- spectral method: 

Find a mapping u N : [0,T] S N such that 

(111. 6) V t, 0 < t < T .Vj.O^j^N , [0u N /dt + (1/2) a N (uJ) + a a 3 u M /ax 3 ](Cj) = 0 , 
Vj,0< j <N , u N (0,Cj) = u°(Cj) • 

or , equivalently, since these equalities entail equalities between polynomials of S N 

(111. 7) [vv€S n , Vt,0<t<T , (du N /at + (1/2) a N (uJ) +aa 3 u N /ax 3 ,v) N = 0 , 

u N (.,0) = I N u° . 

The proof of the existence of the solution to this scheme and of the convergence of u N to u 
will be more technical than the one for the Galerkin method. 

Let us introduce two extended problems that will be useful in the analysis of (II 1.6). The 
first one is a standard K.d.V. problem with initial condition v° ( that will be equal to u(.,t) for 
various times t ): 

av/at + v av/ax + a a 3 v/ax 3 = o , x e ir , t>o , 

(111. 8) v(x+2n,t) = v(x,t) , x € IR , t>0 , 
v(x,0) = v°(x) , x € IR . 

The second problem is a collocation pseudo-spectral problem whose initial condition v^ € S N is a 
suitable approximation of v° , which may differ from i N v° (the one that was taken in (1 1 1.7)): 

Find a mapping v N : [0,T] -*■ S N such that 

Vt,0<UT ,Vj,0<sj«:N , [3v N /at + ( 1 /2) a N (v*) + a a 3 v N /ax 3 ](tp = 0 , 
Vj,0^j<N , v N (O.Cj) = vJ(Cj) , 
or equivalently 

(1 1 1.9) v f € S N , v t, o ^ t < t , (av N /at + (i/2)a N (v„) + a a 3 v N /ax 3 , h/) n = o , 

v N (. ,0) = V N 

We first exhibit a time interval [0, t*] ( tj < T) in which there exists a solution to the 
collocation problem (1 11.9). Then we prove that this solution, together with its three first 
derivatives, can be bounded in [O.t*] by some constants depending only on the initial data. This is 
accomplished in Lemma 1 1 1.2 and in Lemma III. 3. We then prove an estimate of the error between 
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v - v N in terms of v° - v£ . Finally, in section III. 4, we show by an iteration argument that the 
above convergence result, applied to u - u N , can be in fact extended to cover the whole time 
interval [0,T]. 

111.2 Three lemmas about the boundedness of the solution of the collocation 
problem^ 

We begin by noting that, due to the classical theory of differential systems, problem (1 1 1.9) 
admits a local solution. This means that there exists t 0 > 0, such that for all t ^ t 0 , the solution of 

(1 1 1.9) exists and is unique ( note that t 0 may depend on N ). 

However, no information about the boundedness of the solution in any norm independently of 
N is provided from this result. For this we shall take now well suited test functions in (1 11.9) as 
we did in the previous section. 

The first choice ^ = 1 in (III. 9) yields 

(111.10) (a/at) [jo"y N (x,t) dx] = o , 
as in (11.2) , or again 

(III. 1 1 ) V t, 0 < t < t 0 , jo" v N (x,t) dx = Jq" v n (x, 0) dx = jg" v°(x) dx , 
which expresses the conservation of the Fourier coefficient (v N ) 0 (t) . 

We turn now to the proof of 

Lemma 111,1 : For any real number R, there exist three positive constants tj ^ t 0 , ($ 0 , 

depending only on R such that for any initial value v|J verifying 

(IM.12) || v° ||, ^ R 

and any t , 0 < t < tj , we have 

(11113) ||v N (.,t)|Ue 0 (t; -tr 1/2 , 

(111.14) II V„(..t) II. £ Y.(1 + (t* -t)- 5/6 ) 

n " 1 ' * v v 

Proof : Let us take 4 » = v N in (1 11.9); this gives (we drop again the dependence on x and t ) 

( i /2) (a/at)(v N , v N ) N + ( i /2) ((a/ax) i N (v 2 ) , v N ) N + a (a 3 v N /ax 3 ,v N ) N = o . 

From property (1 1 1.4) the third term vanishes by integration by parts, thus 
(a/at) || v N || 2 + ((a/ax) i N (vJ),v H ) = o . 

Using now the Cauchy-Schwarz inequality gives 


i 


(a/at) || v N || 2 ^ |( i N (v 2 ), av N /ax) | < || i N (v 2 ) || ||ev N /ax || 

^ (l N < V N> ’ ^n( V N^N II V N 111 ^ ( y N ’ V N II V N Hi ^ H V N Hl°° H V N H II V N Hi 
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Then by the following Gagliardo-Nirenberg inequality 

(111.15) Vf€H 1 (0,2n) , || tp || L oo < C, || 4 ) || 1/2 !| tp ||; /2 , 
we have 

(111.16) (3/Qt) || v N || 2 < C,( || v N || || v N ||,) 3/2 . 

Using now 41 = l N (v^ + 2a9 2 v N /9x 2 ) in (111.9) yields 

(dv N /9t, v^ + 2a9 2 v N /9x 2 ) N + ( 1 /2)(9 y/9x,\j/) n = 0 . 

The second term vanishes using (1 1 1.4) and integration by parts, hence 

(111.17) (9/9t)[ <x||9v n /9x|| 2 -(1/3)(vJ , 1 ) N ] = 0 . 

Let us integrate this equality between 0 and t, with t < t 0 and use (1 1 1.4) to get 
|| 0v N (. ,t)/3x || 2 - 1| 3v N (. ,0)/9x || 2 = ( 1 /3a)[(vJ(. ,t) , 1 ) N - (v 3 (. ,0) , 1 ) N ] 

^ ( 1 /3oc)(|| v N (.,t) || L oo||v n (. , t) || 2 + ||v N (. ,0) || L oo || v N (. ,0) || 2 ). 
Once more, using again inequality (lll.l 5) gives the following result, 

II av H (. ,t )/ax f < k° * c 2 ( || v N |||" || v N | 5/2 ) , 

where C 2 = ( 1 /3a) C, and 

K° = || 9v n (.,0)/9x || 2 + C 2 || v N (.,0) || L oo || v N (. ,0) || 2 < R 2 ( 1 + C 2 R) . 

Now a bound for the H^(0,2n)-norm of v N (. ,t) is easily recovered by means of (1.9). We deduce 

(111.18) |v N (.,t)||?<K‘ + C 2 (||v N (.,t)|| 1 ,/2 ||v N (.,t)|| 5/2 ) , 
where K' = 2 K° + (J| n v N (x,0)dx) 2 < 2 R 2 ( 1 + C 2 R) + 2n R 2 . 

For some technical reasons which will be clarified at the end of our proof, if necessary we 
take in (III. 1 8) a possibly larger K' in order to satisfy the following inequality 

(111.19) K* > l| v° || 3 . 

This is achieved for instance for K' = max {R 3 , 2 R 2 (tt + 1 + C, R)}. 

Let us sett; = min(T ,1 /(( 1 +C 2 )C 2 K ,2/3 ) and 0 2 = 1 /[( 1 +C 2 )C 2 ] . We show that (111.13) holds 
with these constants. 

Assume by contradiction that there exists t* < t; such that (111.13) is not verified for 
t = f ,i.e. 

(111.20) || v N (.,t*) || 2 > [0/K ') 2/3 - d + c 2 )c 2 t*]- 1 . 

We shall prove now that the derivative of the mapping t -► || v N (.,t) || is bounded, then we will 
derive a lower bound for || v^ || incompatible with (111.19). First, we notice that, for any time 
1 5% t 0 such that 

(111.21) ||v N (.,t) || ^ K' 1/3 , 
we get 
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K- < || v K (. ,t) || J « II v M ( ,t) n; /2 || v N < ,t)|| 5/2 , 

since 

II v M ( .0 II, >||v N (.,t)|| . 

From (111.18) It follows that for any t satisfying (III 2 1), we have 

(111.22) ||v N (.,t)||^(HC 2 )(||v N (.,t)|| i 1/2 ||v N (. I t)|| 5/2 ) ; 
we deduce then 

(111.23) || v N (.,t) ||, < ( 1 +C 2 ) 2/3 || v N (.,t) ll 5/3 
Noting that 

[( 1 /K’) 2/3 - ( 1 +C 2 )C 2 t*]” 1 > k ,2/3 , 

we deduce from (111.20) that (111.21) holds for t = t* , hence (111.23) is satisfied for t = t*. 

Let us introduce now (111.23) in (1 1 1.1 6). For any t, 0 < t < t* such that || v N (.,t) || ^ K' 1/3 we 
have 

(a/at) II v N (.,t) II 2 < ( 1 + C 2 )C 2 II v N (.,t) || 4 , 
which can be written as 
(IM.24) (a/at) [-i/||v n || 2 ]^(uc 2 )c 2 
Let us consider now the set Ac [ 0 ,t* ] defined by 

A = {s, 0 < s < t* such that for any t in [s ,t* ] we have || v N (. ,t) || ^ K' 1/3 } . 

It is an easy matter to check that, by virtue of the continuity of the function t -*■ || v N (. ,t) || , A is a 
closed interval [a*,t*] of [0,t*]. Besides, from (111.24) we obtain for any s in A that 
II v N (.,s) || 2 £ [( 1 /|| v N (. ,t») || 2 ) ♦ ( 1 + C 2 )C 2 (f-s)]' 1 , 
and by (111.20) 

(111.25) V s € A , ||v N (.,s)|| 2 >[(1/K , ) 2/3 -(UC 2 )C 2 s]- 1 > K‘ 2/3 . 

Applying (111.25) to a* shows that o* = 0 hence A turns out to be equal to [0,t*]. We arrive now 
to the contradiction between (111.25) with s = 0 and (111.19). We deduce first that the max time t. 
is independant of N since no explosion occurs before tj and that (111.20) cannot hold. Whence, we 
have (111. 13) for all t < tj!| .Now (111.14) follows from (111.13), (111.18) and Cauchy-Schwarz 
inequality. 

Remark III. 2 : We note that equation (111.17) is the discrete analogous of the equation (11.4) of 
conservation of the third energy integral for the solution of the collocation problem. 
Unfortunately, the conservation of the second integral does not hold any more in the current case. 
This failure does not allow us to deduce directly the conservation of the H^(0,2n)-norm of v N . On 
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the other hand, this would still be possible if the following skew symmetric decomposition of the 
nonlinear term 

(1/3) 0 N (v 2 ) + (1/3) v M 3v N /3x 

was considered in (1 1 1.9) instead of (1/2) 3 N (vj|j). However, this choice would slow down the 
efficiency of the numerical scheme since an extra nonlinear term should be computed. Besides, we 
will prove also for the genuine pseudospectral scheme (1 1 1.9) a uniform bound of the 
H^(0,2n)-norm of v N on the large time interval [ 0 ,T ] . To this end, we start by proving a bound 
on a "small" time interval in the next corollary. 

Corollaru 111.1 : For any real number R, there exist three positive constants t, , 
depending only on R such that for any initial value vjj verifying 

(111.26) || vj H, <£ R 

and any t , 0 < t < t , , we have 

(111.27) || v„(.,t) || < e,, 

(111.28) ||v N (.,t) II, < . 

Proof : It is a simple consequence of the previous lemma. For instance one can take 1 1 = tj /2, in 
which case (111.27) and (111.28) hold with (5, = § 0 [2/to] 1/2 and =y Q ( 1 +( tJ/2)" 5/6 ) 

The same kind of bounds we obtained for v N will now be proved for the derivative 3v N /dt. Let 
us first recall that in view of (111. 1 0) the following equality holds 

(111.29) Vt,0<Ut 0 , fa" 0v N /at(x,t) dx= 0 . 

Further information concerning the boundedness of the L 2 (0 ,2 tt) and H^(0,2n)-norms of 3v N /3t 
are obtained by differentiating equation (II 1.9) with respect to the time variable. We get 
(in.30) v y € s n , v t, o <t <t , o 2 v N /at 2 + a N (v N av N /at) + a (a 3 /0x 3 )(0v N /at), <v) N = o . 

We can now prove the following result 

Lemma III. 2 : For any real number R, there exist three positive constants t* ^ t, , and ** 
depending only on R such that for any initial value v£ verifying 

(111.31) ||vJ|| 4 <R 

and any t , 0 < t < t* , we have 

(111.32) || (0v N /3t)(.,t) || 6*. 

(111.33) H (av N /3t)(..t) ||, . 
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Proof : Let us take y □ 9v N /9t in (111.30). We obtain, using (III. 4) and noting that 9v N /9t is 
periodic 

V t, 0 < t < t, , (1/2)(9/9t)||9v N /9t|| 2 + (9 N (v N 9v N /9t),9v N /9t) N = 0 . 

By the definition (III. 2) and the property (III. 4), integrating by parts and using the 
Cauchy-Schwarz inequality, it follows 

V t, 0 < t < t, , (9/9t) || 9v N /9t || 2 < 2 |( I N (v N 9v N /9t),(9/9x)9v N /9t)| 

<2||l N (v N 9v N /9t)||||9v N /9t|| 1 , 

< 2 II v N |t L oo 11 9v N /9t 11 II 9v N /at ll, . 

Using the Oagliardo-Nirenberg inequality (III. 15) and Corollary III. 1 , we can find a constant K" 
depending only on the H^(0 ,2n)-norm of the initial condition such that 

(111.34) Vt.O^Ut, , || v N || L oo ^ K" . 

We can therefore obtain the inequality 

(111.35) Vt,0<Ut, , (9/9t) || 9v N /9t || 2 < 2 K" || 9v N /9t || || 9v N /9t ||, . 

Our goal is now to provide a bound for || 0v N /9t ||, . To this end, let us take 

4 t = I N (v N 9v N /9t + a(9 2 /9x 2 )9v N /9t) 
as test function in (111.30). This choice yields the equality 

-(9 2 v N /9t 2 ,a(9 2 /9x 2 )9v N /9t) = (9 2 v N /9t 2 , v N 9v N /9t) N . 

After integration between 0 and t, 0 < t ^ t v we obtain 

(111.36) || (9/9x)(0v N /9t)(.,t) II 2 - II (9/9x)(9v N /9t)(.,0) || 2 = 2a _1 Jo (9 2 v N /9t 2 , v N 9v N /9t) N ds. 
Let us now focus on the right-hand side of the previous equality. Integrating by parts with respect 
to the time variable gives 

lo(9 2 v N / 9t 2 , v N 9v N /9t) N ds n [(9v N /9t, v N 9v N /9t) N ]q - jo (9v N /9t,(9/9t)(v N 9v N /9t)) N ds 

= [(0v N /9t, v N 9v N /9t) N ]J - Jo (9v N /9t,(9v N /9t) 2 ) N ds 

_ f l (Ply, v a 2 v/ /at 2, l He 

JO ' N' ” ' 1 ’ N “ ’N' • 

Noting that the last term on the right-hand side is the opposite of the left-hand term , we get 
2 Jo(9 2 v N /9t 2 , v N 9v N /9t) N ds = ((9v N /9t) 2 (. ,t) , v N (. ,t)) N - ((9v N /9t) 2 (. ,0) , v N (. ,0)) N 

_ Jo( 1 X9v N /9t) 3 ) N ds 

By virtue of (111.34) and the Oagliardo-Nirenberg inequality (111.15), we deduce that 

2 IJo0 2 V at2 , v n av N /at ) N ds k K" || (9v N /9t)(. ,t) || 2 + II (9v N /9t)(.,0) || 2 || v° || L « 

+ Jo II 9v N /9t || L oo || 9v N /9t || 2 ds , 

^ K" || (9v N /9t)(.,t) || 2 + || (9v N /9t)(.,0) || 2 || vj || L oo 
+ C, Jo II av N /9t ||j /2 II 9v N /9t || 5/2 ds . 
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Injecting this inequality in (111.36) and using (1.9) and (111.29) gives 

(a/2) || (0v N /0t)(. ,t) II? ^ a || (0v N /Ot)(.,O) ||? + || (0v N /0t)(. .0) || 2 || vj | L « 

+ c, ll II av N /at||, ,/2 1| av N /at || 5/2 ds + k- || (av N /at)(.,t) || 2 . 

We can rewrite (1 1 1.9) as follows 

av N /at = ( i /2) a N (v 2 ) + a a 3 v N /ax 3 . 

Then using (1 1 1.3 1 ), it is an easy matter to deduce from this equality on the one hand that 
||(0v N /0t)(.,o)H c [ ||a/ox((v2) - i N (v 2 )) || - ||a/ax(v 2 )) || + 1| a\/ax || 
and using (111.5) with r = s = 1 , we deduce that ||(av N /at)(. ,0)|| is bounded by R 2 ( 1 + a R) ; on the 
other hand, by taking its derivative in the x-direction, we derive that the term || (8v N /at)(. ,0) fl, 

O 

can be bounded by R (1 + a R). Finally , with a new constant K" depending only on R , we derive 
the inequality 

(111.37) Vt.CUUt, , ||(8v N /at)(. > t)||?<K"(1 +l|(0v N /at)(.,t)» 2 ) 

+ cj‘ l|av N /at|i; /2 ||av N /at|| 5/2 ds . 

Let us now consider the set B defined by 

B = { s € ] 0 , t 1 ] : V t € [0 ,s] : C, j‘ 0 || 9v N /3t ll" 2 || 3v N /at || 5/2 ds < K" ) . 

It is an easy matter to check that B is not empty; indeed, we derive from (1.8) and (1 11.35) that 

, (a/at)||av N /at(.,T)|| 2 < 2cNR"||av N /at(.,x)|| 2 . 

Then, we get, using again (1.8) 

v x, o ^ x < f \ || av N /at(.,x)|| 2 + cn -1 || av N /at(.,x)||? < 2 1| av N /at(.,o) || 2 exp(2cN k m x). 
This proves the existence of a time t { > 0 (obviously depending on N) such that ]0 ,tj] c B. Next, we 
derive from (111.37) that 

(in.38) v t € b , ||(av N /at)(.,t)||?< K" (2 + ||(av N /at)(.,t)|| 2 ) , 
and due to (111.35) 

v t € b , (a/at) || av N /at || 2 «sK"yic( 2 1| av N /at || + || av N /at || 2 ) 

<4K"v r K 7r (u||av N /at|| 2 ) . 

After integration on time, we obtain 

(111.39) V t € B , (2 + ||(3v N /3t)(. ,t)|| 2 ) < ( 2 + ||(8v N /3t)(. ,0)|| 2 ) exp(4 K" -/K 77 ” t) 

From (111.38) we then state 

(111.40) V t € B , ||(8v N /at)(.,t)||? < K" ( 2 + ||(av N /at)(.,0)|| 2 ) exp(4 K" yrt) 

We deduce that there exists a constant x 0 , that depends only on the H^(0,2n)-norm of the initial 
function v^ and on T such that 

Vt€B , ||(av N /at)(.,t)||?<x 0 . 
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From the definition of B we can state that this set contains an interval [ 0 , t* ] with 
t* = min( t , , K’VCCji^ 2 ), which is therefore independent of N. Thus the desired results follow 
from (111.39) and (111.40). 

In order to prove the convergence of the discretization we need a further stability result 
l emma III. 3 : For any real number R, there exists a constant * 3 depending only on R such that 
for any initial value v£ verifying 

(111.41) ||v°|| 4 <R 

and any t , 0 < t ^ t* , we have 

(111.42) || v M (.,t) || 3 < * 3 . 

Proof : Let us choose y = 0 3 v N /9x 3 in the equation (1 1 1.9). We obtain 

| oc 1 1| 9 3 v N /0x 3 || 2 = -(0v N /0t,0 3 v N /9x 3 ) - ( 1 /2)(0/9x( I N (vJ)),a 3 v N /0x 3 ) 

^ II av N /at || ||a 3 v N /ax 3 1| + [|| a/ax((vj)- i N (vJ)) || + c || v N || 2 ] || a 3 v N /ax 3 1| . 
Using ( 1 1 1 .5) with r □ s = I we have 

(Ml 43) | a | II a 3 v N /ax 3 II 2 < II av N /at II ||a 3 v N /ax 3 II + C ||v N ||* ||a 3 v N /ax 3 || . 

On the other hand, using the definition of the H 2 (0,2n)-norm , the inequality (111.15) and the 
inverse inequality (1.8) we deduce 

b v* i a < o n v M n 2 h v H a, * n v B g?] <cniivX 

By virtue of (111.43) we obtain 

lex in a\/jx 3 || J < || av„/at || ||a 3 v H /ax 3 1| * c || v„ \\] || a 3 v„/ax 3 1| . 

The two previous lemmas yield now that the term || o v H /3x 3 1| is bounded and we derive the desired 
result (1 11.42) by noting that 

II a 2 v N /ax 2 Nc|| a 3 v N /ax 3 ||. 

111.5 Analysis of the convergence in the local interval [0,t*]. 

We are now in the position of stating a first local convergence result of the solution of the 
collocation problem (III. 9) to the one of the K.d.V. equation in the time interval [0,t*]. Precisely 
we have 

Proposition 111.1 : Assume that v° belongs to H™(0,2n) , for some m ^ 4 , and that vj] is 
bounded in (0,2n) independently of N. Then there exists a constant A m > 0 , depending 
continuously on || v£ || m and independent of N such that for any t, 0 < t < t*: 

(111.44) | v N <„t) - vC..t) H, <A m N 2 -"'-K-v°||,. 
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Proof : From (1 1 1.8) it is easy to check that (remind that E(x) denotes the integral part of x) 

V t , 0 < t < T , P E(N/2) (8v/at + v 8v/3x + a a 3 v/3x 3 ) = 0 , 

so that 

V y € L 2 (0 ,2n) , V t , 0 ^ t<T , (P E(N/2) (8v/8t + v 3v/8x + oc 3 3 v/8x 3 ) .y) = 0 . 

For a fixed time t, let us substract this equation from the collocation equation (1 1 1.9) at the same 
time. If we set v(t) = P E(N/2) v(t) and e(t) □ v(t) - v N (t) , for any q# in S N we derive the identity 

(111.45) (9e/at + a8 3 e/ax 3 ,^) = (1/2) (a/3x(l N (v 2 ) - P E(N/2) (v 2 )), <jr) . 

Taking \j/ = e(t) for all t and noting that v 2 (t) € S N we obtain 

(111.46) (a/at) || e|| 2 = (8/8x(v 2 - v 2 ) ,e) + (8/8x(v 2 - P E(N/2) (v 2 )) ,e) +(a/8x(l N (v 2 - v 2 )),e). 
The last term on the right can be bounded as follows 

(8/8x(l N (v 2 - v 2 )),e) 3 ((v 2 - v 2 ),8e/8x) N 3 ((2v - e)e,8e/ax) N 

2 (v e,v e)" 2 II 8e/8x || + || 3e/3x || L oo || e || 2 
< 2 || v || L oo || 3e/3x || || e || + || 8e/3x || L oo || e || 2 . 

Then from (111.46), we deduce 

(1 1 1.47) (3/at) || e || 2 < || e ||{[|| 8v/8x || L oo + || 3v/8x || L « ] || v - v ||, + || v 2 - P E(N/2) v 2 », 

+ 2 || v \\ L oo || 8e/8x || + || 3e/3x || L oo || e ||} . 

From (1.2) and (1.7) we get, as soon as v° belongs to H™(0,2n) with m ^ 3 

(111.48) [|| 3v/8x || L cc + || 8v/ax|| L «o ] || v - v ||, +|| v 2 - P E(N/2) v 2 ||, < C(m) N 1 '" 1 , 

Here C(m) is a constant that depends only on the norm i} m || v° || m . Moreover, from (111.42) we 
have (we enlarge, if needed, the value of C(m) keeping although its dependence only on T) m || v° || m ) 

(111 .49) 2 || v || L oo + || ae/8x || LOO < C(m) + ^3 . 

Taking into account (111. 47) (111.48) and (111.49) 

(111.50) (3/at) || e || 2 < (C(m) + * 3 ) (N 1_m || e|| + || 8e/3x || || e|| + || e|| 2 ) , 
where the constant depends on the initial conditions only. 

Our next goal is to obtain an estimate for || 8e/ax || in order to use it in (111.50). For this, 
we take iv 3 (i N vjjj - P E ( N/2 )V 2 ) - 2aa 2 e/3x 2 in (III. 45) , so that we find the identity 
-2oc(ae/3t,3 2 e/8x 2 ) = (8e/at,P E(N/2) v 2 - I N vJ) . 

Therefore 

(111.51) a(8/3t) || 8e/8x || 2 3 (3e/at, P E(N/2) v 2 - v 2 ) + (3e/at, I N (v 2 -vJ)) , 
or again , if we set w 3 P E(N/2) v 2 - v 2 , 

(111.52) oc(a/at) || 3e/3x || 2 3 (3e/at,w) + (de/at, l N [e(v + v N )]) . 


Besides we remark that 
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(3e/3t, I N [e(v + v N )]) = (3e/0t,e(v + v N )) N = ( 1 /2)(3e 2 /0t,v + v N ) N . 

Then integrating (111.52) between 0 and t for any t < tj we derive 

(111.53) oc|| (ae/axX.,0 || 2 - cx|| (ae/ax)(.,0) || 2 = Jq ( 3e/at,w) ds +(1/2)jJ, (3e 2 /3t,v + v N ) N ds . 
Setting z = v + v N and integrating by parts with respect to t gives 

(111.54) | Jo(0e 2 /3t,z) N ds |= | -[(e 2 ,z) N ]^ + (e 2 ,az/at) N | 

< | (e 2 (.,t),z(.,t)) N l + |(e 2 (.,0),z(.,0)) N | + max [ ||3z/at(. ,x)|| L oo ]}q || e|| 2 ds 

0<x<t 

Furthermore, we have 

(111.55) | fo(3e/at,w) ds | = | -(e(t),w(t)) + (e(O).w(O)) + ^(e.aw/at) | 

< II e(t) HU w(t) || + || e(0) IHI w(0) || + t max || e(x) II, || 3w/3t(x) ||_, 

0 < x < t 

Finally, we can deal in a similar way the last term of the right-hand side of (III .53). 

Let us set now 

K 0 r || 0e/0x(O) || 2 + | (e 2 (0) ,z(0)) N | + || e(0) || || w(0) || ; 

then, from (III. 53) to (111.55) we deduce that t 

V t, 0 < t < t* ,| 3e/3x(t) || 2 K 0 + c max {[ || z(x) || L <*> + || 3z/3t(x) || L ~] J 0 || e(s) || 2 ds 

0 < x < t 

+ || 3w/3t(x) ||f, + (1/2c) || e(x) || 2 } 

+ II e(t) || 2 + || w(t) || 2 . 

It is an easy matter to derive from (1.2) and (III. 8) that if v° belongs to H™(0,2 ti) with m > 3, 
then v belongs to L“(0,T; H™(0,2n)) and 3v/at belongs to L“(0,T; H™ _3 (0,2 tO). Thus, by a 
straightforward application of (1.7) , we deduce 

K 0 + C max [|| 0w/at(x) || 2 , + || w(x) || 2 ]< C(m) {N 2(2_m) + || e(0) || 2 } . 

0 < x < t* 

Hence we derive from (111.28) and (111.33) that 

(111.56) max || 0e/0x(x) || 2 < C(m)[N 2(2_m) + J 0 || e(s) || 2 ds +||e(0)|| 2 ] 

0 < x s* t 

Let us now introduce this last inequality into (111.50); we obtain 

(3/at) || e(t) || 2 (C(m) + * 3 ) (N 2<2 ' m) + || e(t) || 2 + {* 1| e(s) || 2 ds + || e(0)|| 2 ) , 

whence, after integration with respect to the time variable 

Vt,0<Ut; ,|| e(t) II 2 < II e(0) || 2 + (C(rn) + * 3 )(t(N 2(2 - m) + || e(0) || 2 ) + f 0 || e(s) || 2 ds) 

< (C(m) +'jf 3 )(t(N 2<2-m) + || e(0)|| 2 ) + JqII e(s) || 2 ds ) . 


The Gronwall lemma yields now 
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dll. 57) Vt,0<ut; , ||e(t) II 2 < (C(m) + Y 3 )(N 2<2 ‘ m, + ||e(0)||2 ) 

We recall that the constant (C(m)+Tf 3 ) that appears here depends only on the H™(0 ,2n)-norm of 
the initial condition v° and on *3 , which in turn, from Lemma 1 1 1.3, is a bound for the 
H^(0,2n)-norm of v N . Hence c is a constant independent of N and t.. 

The result (111.44) follows now using (1.7) once more together with (111.56) and (111.57). 


111.4 Convergence results for the approximation 


We can state now the main result of this section. 

Theorem III. 1 Assume that u° belongs to H™ (0,2rt) , for some m > 4 . Then for anyt, 0 < t < T 
and any N large enough the following estimate holds 
(111.58) ||u(.,t)-u N (.,t) ||, < c N 2-m 

Proof: Let us choose a class of initial conditions for ( 1 1 1 .9) such that 
(11159) II vj || 4 < 2 tj 4 II u° || 4 , 

where *g 4 is the constant which appears in (1.2). From Lemma II 1.2 ( applied with 
R = 2 || u° || 4 ) , we deduce the existence of a time t* independent of N such that, for any initial 

condition verifying (111.59), the solution of (1 1 1.9) exists for any t , 0 < t < t* . Under the 
current hypotheses, by virtue of Proposition III. 1 there exists a constant A m = A( || u° || m ) such 
that, for any time t , 0 < t < t*. 

(111.60) ||v N (.,t)-v(.,t)|| 1 <A m N 2 - m + ||v;-v°||, . 

Finally, we denote by N* the integral part of 1 + ([(4K + 1 3)A m ] -1 i] 4 || u° || 4 ) 1/<s " m> , where we 
have set K = T/t*. 

From (1 1 1.5), it can be deduced that there exists a constant N, ^ N* such that for any 
N ^ N, the estimate (111.59) is true for u^ and that 


(111.61) ||u°-uj||, + ||u(.,t)-P N u(.,t)||, + N" 3 || u(. ,t) - P N u(. ,t) || 4 < A m N 1-m 
(if necessary, we have enlarged the value of the constant A m ). 

We are going to prove by induction on k ^ K. that 

For any t , kt* < t ^ (k+ 1 )t* , u N exists and satisfies 
(H k ) II u N (. ,kt*) || 4 < 2 tq 4 || u° || 4 and 

V t ,kt* <t<(k+1)t* , || u(.,t) - u N (.,t) ||, ^ (k + 2)A m N 2 ‘ m 
It is an easy matter to check that (H 0 ) is a simple consequence of (111.60) and (111.61 ). Let us 
assume that H k is true for k and let us prove it for k+1. It is readily seen that the solution 
of ( 1 1 1 .8) with v° = u(. ,(k+ 1 )t*) is the solution of (1.1) and that the solution of (1 1 1 .9) with 


- 22 - 


v n = u N (. ,(k + 1 )t*) is the solution of (111.7) for t > (k + 1 )t*. First, we have to prove a bound for 
||u N (. ,(k+ 1 )t*)|| 4 . Using (1.7) and the inverse inequality (1.8) we deduce from (111.61) and the 
previous estimate that 

Vt,0<t<(k + 1)t* , ||u N (.,t)-P N u(.,t)|| 4 ^4N 3 ||u N (. 1 t)-P N u(. 1 t)|| 1 

<4N 3 [|| u N (.,t) - u(,t) ||, + ||u(.,t)-P H (.,t) ||,] 
< 4 (k + 3)A m N 5-m . 

Therefore 


V t ,0 <U(k + 1)t* , ||u N (. I t)|| 4 <||u N (. > t)-P N u(.,t)|| 4 + ||u(. I t)-P N (. I t)|| 4 + ||u(.,t)|L 


^(4K + 1 3) A m N 5_m + || u(. ,t) || 4 . 

If N is chosen greater than N* we derive from (1.2) that || u N (. ,(k + 1 )t*) || 4 < 2 *g 4 || u° || 4 
Moreover, from Lemma 111.1 , we deduce that u N exists for any t , (k + 1 )t* < t < (k+2)t* and 
from (111.60) 


Vt,0<Ut* , || u N (. ,(k+ 1 )t*+t) - u(. ,(k+ 1 )t*+t) ||, < A m N 2_m 

+ II u N (.,(k + 1 )t*) - u(.,(k+ 1 )t*) ||, . 

Using now the induction hypothesis we obtain 

vt,0<ut* , || u M (.,(k + 1)t*+t) - u(.,(k+1)f+t) II, < (k + 3)A m N 2 - m . 

This proves hypothesis (H k+1 ). Thus the induction procedure applies successfully and the desired 
result (111.58) follows. 


Remark III. 3 : The stability result in the H 4 (0,2n)-norm stated in (H k ) plays a fundamental 
role in the proof of the global convergence result (111.58). As we have seen, stability has been 
obtained by exploiting the spectral decay of the error on each local interval [kt*, (k + 1 )t*]. 
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